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ABSTRACT 

Context. More and more evidence lead to consider classical Be stars as stars rotating close to the critical velocity. If so, then the 
question which naturally arises is the origin of this high surface velocity. 

Aims. We determine which are the mechanisms accelerating the surface of single stars during the Main Sequence evolution. We study 
their dependence on the metallicity and derive the frequency of stars with different surface velocities in clusters of various ages and 
metallicities. 

Methods. We have computed 112 stellar models of four different initial masses between 3 and 60 Mq, at four different metallicities 
between and 0.020, and with seven different values of the ratio fi/f\. rit between 0.1 and 0.99. For all the models, computations were 
performed until either the end of the Main Sequence evolution or the reaching of the critical limit. 

Results. The evolution of surface velocities during the Main Sequence lifetime results from an interplay between meridional circu- 
lation (bringing angular momentum to the surface) and mass loss by stellar winds (removing it). The dependence on metallicity of 
these two mechanisms plays a key role in determining for each metallicity, a limiting range of initial masses (spectral types) for stars 
able to reach or at least approach the critical limit. Present models predict a higher frequency of fast rotating stars in clusters with 
ages between 10 and 25 Myr. This is the range of ages where most of Be stars are observed. To reproduce the observed frequencies 
of Be stars, it is necessary to assume first that the Be star phenomenon occurs already for stars with v/v alt > 0.7 and second, that 
the fraction of fast rotators on the Zero Age Main Sequence is higher at lower metallicities. Depending on the stage at which the star 
becomes a Be star, the star at this stage may present more or less important enrichments in nitrogen at the surface. 
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1. Introduction 

The aim of the present work is to study how the surface veloc- 
ity varies as a function of time during the Main Sequence (MS) 
evolution for stars of various initial masses, metallicities and ro- 
tation. Although predictions of the evolution of the surface ve- 
locity can already be found in literature (see Heger & Langer 
2000; Meynet & Maeder 2000, 2005), the published results lack 
homogeneity in their physical ingredients and cover too narrow 
ranges of initial velocities to provide valuable answers to the fol- 
lowing questions: 

- When is very fast rotation encountered in the course of the 
evolution of massive single stars on the MS phase? 

- What is the percentage of stars rotating faster than a given 
value of v e q/v cr i t in stellar clusters of various ages and metal- 
licities? 

- What are the expected chemical enrichments at the surface 
of these very fast rotating stars? 

- What are the initial conditions required for single stars to 
reach the critical velocity limit (defined as the rotational ve- 
locity such that the centrifugal acceleration exactly balances 
the gravity at the equator)? 

- Do these conditions depend on the mass and metallicity and 
if yes how? 

To know which are the stars which can encounter the critical 
limit during the MS evolution is interesting in at least two astro- 
physical contexts: first, the study of Be stars, second, the very 
metal-poor stars representative of the stellar generations in the 
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early Universe. Let us briefly describe the link with these two 
research areas. 



1.1. Be stars and the critical limit 

It is now widely accepted that the "Be phenomenon" is linked 
to fast rotation. Indeed, classical Be stars are B-type stars that 
exhibit line emission over the photospheric spectrum. This emis- 
sion probably arises from an outwarding equatorial disk (Porter 
& Rivinius 2003) whose formation is (at least in part) due to the 
fast rotation of the star. If the star is rotating with v e q/fcrit ~ 
0.95, the additional velocity necessary to launch material into 
orbit is of the same order as the speed of sound in the outer 
layers (Townsend et al. 2004). Therefore many instabilities, as 
e.g. non-radial pulsations, begin to be effective for orbital ejec- 
tion (Owocki 2004). Very interestingly, recent interferometric 
observations brought strong support to the idea that these stars 
are very fast rotators: the oblateness of the Be star Achernar 
(Domiciano de Souza et al. 2003) suggests indeed that the star 
is rotating near the critical limit; Meilland et al. (2007) obtained 
for the first time direct evidence that the disk around the Be star 
a Arae is in Keplerian rotation. The wind geometry is compat- 
ible with a thin disk wind and a polar enhanced wind. These 
authors also found that a Arae is rotating very close to its criti- 
cal rotation. Even the mystery behind the ji Ceph case seems to 
be elucidated: /3 Ceph is a known Be star but also a known slow 
rotator, thus challenging the fast rotation scenario. However f3 
Ceph is a multiple system. Schnerr et al. (2006) have shown that 
/? Ceph, the slow rotator, produces most of the light but is not 
responsible for the emission lines. The emission lines originate 



2 



Ekstrom et al.: Origin of Be stars 



from a much fainter companion, the Be star, which they consider 
as probably fast rotating. 

If the above observations do support the view that Be stars 
are at the critical limit, there are still some doubts that all Be stars 
are indeed at this limit. For instance Cranmer (2005), on the basis 
of the observed distribution of the v sin i for Be stars, conclude 
that while late type Be stars could be at the critical limit, early Be 
type stars might rotate at only 40% of their critical velocity. Let 
us note that the interpretation of the observations is stained with 
some theoretical uncertainties. A good example is the way the 
deduced velocities are affected by the fact that for a fast rotating 
star most of the radiative flux originates from the polar regions 
(von Zeipel 1924). The low velocity regions, near the pole, will 
thus have a stronger impact on the spectrum than the dimmer 
faster regions near the equator. This implies that the rotational 
velocity of very fast rotators have been probably systematically 
underestimated (Collins 2004; Townsend et al. 2004). Cranmer 
(2005) has accounted for this effect in his work, but the problem 
is sufficiently complex that the results are uncertain. 

Anyway, at the critical limit or not, Be stars represent won- 
derful laboratories for studying the effects of extreme rotation. 

1.2. Very metal poor stars and the critical limit 

The link between metal-free or very metal-poor stars and the 
question of the evolution of the surface velocity certainly de- 
serves great attention, since it might deeply affect the evolution 
of these stars. Let us recall that at low metallicity, the radiation 
driven stellar winds are much weaker than at high metallicity, 
thus it is commonly accepted that mass loss for very metal-poor 
stars is very inefficient. This might not be true however when ro- 
tation is taken into account in the models for at least two reasons 
(see Meynet et al. 2006b, to get a more detailed discussion of 
these points): first and somewhat paradoxically, since the stars 
are losing very little amount of mass by radiatively driven stel- 
lar winds, they have more chance to lose large amount of mass 
by rotational ejection during the MS evolution. Indeed, because 
of the low line driven winds, they will not be able to lose their 
angular momentum. They are thus more prone to reach the criti- 
cal limit during their evolution (Maeder & Meynet 2001 ; Meynet 
et al. 2006a). At the critical limit, as in the Be case, stellar matter 
is launched into an equatorial disk which will eventually dissi- 
pate by radiative effects. In that case the material will be lost by 
the star. If correct, such an evolution would thus promote mass 
loss through a mechanical wind triggered by rotation. This "me- 
chanical wind", as stressed above, appears to be more important 
in metal-poor regions where the radiation driven stellar winds 
are expected to be very low. Second, rotational mixing enhances 
the surface CNO content which will lead to some increase of the 
outer layers opacities. This may revive the line driven winds and 
lead to strong mass loss. 

Such mass loss, either triggered by rotational ejection or 
through the enhancement of the CNO content of the outer lay- 
ers, have (if realized) important consequences on the evolution 
of the first stellar generations as well as on their nucleosynthe- 
sis (Ekstrom et al. 2006; Hirschi 2007; Decressin et al. 2007). 
Also such mass loss may be important in the process leading to 
the formation of a collapsar which is the most favoured model 
for the progenitors of the long Gamma Ray Bursts (Woosley 
1993; Yoon & Langer 2005; Woosley & Heger 2006; Meynet 
& Maeder 2007). 

In the present work, we use rotating stellar models to address 
the questions of: 



- the possible origin of Be stars from a purely theoretical point 
of view: can models reproduce the way the number of Be 
stars varies with the metallicity or with the age? Is their 
number sufficient to account for the observed number of Be 
stars (supposing that Be stars are star rotating with values of 
v eq/^crit superior to a given value)? 

- the behaviour of rotating stars at very low or zero metallicity: 
how does the lack of metals influence the evolution of the 
equatorial velocity? 

Of course the answers given to these questions will depend 
strongly on the stellar models used. However we think that the 
exercise deserves to be undertaken for the following reason: the 
stellar models we use were computed with physical ingredients 
that made predictions consistent with many observational fea- 
tures of massive stars (see Meynet & Maeder 2005, and the ref- 
erences therein on the series of papers by the Geneva group on 
rotating models) and thus it is interesting to explore in more de- 
tail what these models say about fast rotating stars. For very fast 
rotation, the present models probably need still to be improved 
by the incorporation of additional physics (e.g. pulsation), and 
we hope that the results obtained here will guide us towards im- 
provements in the future. 

The paper is organized as follows: in Sect. 2 the physical 
ingredients of the models are presented. We also recall some ba- 
sic relations from the Roche model useful for discussing some 
properties of our models on the ZAMS. In Sect. 3, we present 
various characteristics of our models on the ZAMS as a function 
of the initial mass, velocity and metallicity. The evolution of the 
models on the MS is presented in Sect. 4. The question of the 
initial conditions required to reach a given value of f eq /v C rit dur- 
ing the MS phase is the subject of Sect. 5. The frequency of stars 
with velocities exceeding a given value of Veq/vcnt is presented 
in Sect 6. The sensitivity of the predicted frequencies on various 
hypotheses made in this work are discussed in Sect. 7. Possible 
links with observational features of Be stars will be presented 
in Sect 8. Conclusions and future perspectives are the subject of 
Sect. 9. 

Let us mention that electronic tables are available for all our 
ZAMS and evolutionary models at http://obswww.unige. 
ch/Recherche/evol/Critical- limit- and-Be- stars 

2. Physics of the models 

2.1. Input physics 

Since the detailed description of the physics of the models has 
already been presented in detail in previous papers (see refer- 
ences below), we shall keep here the discussion very short. The 
physics of rotation (transport of the angular momentum and of 
the chemical species) are treated as in Meynet & Maeder (2005). 
These models are based on the theory of the transport mecha- 
nisms induced by rotation proposed by Zahn (1992) and further 
complemented by the works of Talon & Zahn (1997) and Maeder 
& Zahn (1998). 

We have computed the evolution of five different masses, 
M = 1, 3, 9, 20, 60 Mq, at four different metallicities, Z = 0.020, 
0.002, 0.00001, 0. For each combination of M - Z, we have 
computed models with seven different initial rotational rates, 
(o = n/Q crit = 0.10, 0.30, 0.50, 0.70, 0.80, 0.90, 0.99. The 3 
to 60 Mq models have been evolved until either the beginning 
of core helium burning or the reaching of break-up limit. For the 
1 Mq models, only ZAMS models are presented. 

The initial composition of the models for the metallicities 
Z = 0.020, 0.002, 0.00001 and are respectively X = 0.7050, 
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0.7545, 0.7599 and 0.76 and Y = 0.275, 0.2435, 0.24002 and 
0.24. For the heavy elements, we have used the same mixture 
as the one used in the OPAL opacity tables: solar mixture of 
Grevesse & Noels (1993) for Z = 0.020 and 0.002 and a- 
enhanced mixture of Weiss (1995) for Z = 0.00001 The nu- 
clear reaction rates are from the NACRE database 2 . 

The convection is treated according to the Schwarzschild cri- 
terion. The size of the convective core is increased by a moderate 
overshooting of 0.1 H p . 

The 3 Mq models were computed without mass loss. The 9 
Mq models were computed with de Jager et al. (1988) prescrip- 
tion. For the higher mass models we have used the prescription 
proposed by Kudritzki & Puis 2000 (hereafter KP00), and de 
Jager et al. (1988) in the temperature domain not covered by 
the first prescription. The effects of rotation on the mass loss 
rates are taken into account as explained in Maeder & Meynet 
(2000a). 

2.2. The Roche model 

In all the derivations, we shall use the Roche model, i.e. we ap- 
proximate the gravitational potential by GM r /r where M r is the 
mass inside the equipotential with a mean radius r. The radius r 
which labels each equipotential is defined by r = (V r /(4-/3n)) 1/3 
where V r is the volume (deformed by rotation) inside the equipo- 
tential considered. 3 

In the frame of the Roche model and for a conservative rota- 
tion law (i.e. such that the centrifugal force can be derived from a 
potential), the shape of a meridian at the surface of an equipoten- 
tial is given by couples of r and 6 values satisfying the following 
equation: 



GM r 1 „ ., GM r 

+ -fi 2 r 2 sin 2 6 = -, 

r 2 r n 



(1) 



where r is the radius at colatitude 6, fi r the angular velocity at 
that radius, M r the mass inside the considered equipotential sur- 
face and r p , the radius at 6 = 0. Thus the shape of a given equipo- 
tential is determined by three parameters M r , Cl r and r p . The first 
two M r and fi r are independent variables. The third one is de- 
rived from the first two and the equations of stellar structure. 

Setting x = (^r) r, one can write Eq. (1) (see Kippenhahn 
& Thomas 1970) 



- + -x z sin z 6» = — . 

X Is X-p 



(2) 



With this change of variable, the shape of an equipotential is 
uniquely determined by only one parameter x p . 

With M the total mass of the star and fi the angular velocity 
at the surface, Eqs. (1) or (2) give the shape of a meridian line at 
the stellar surface. Setting 



R e 

/= — , 



(3) 



1 see http : //www-phys . llnl . gov/Research/OPAL/opal . html 
for exact references 

2 see http : //pntpm . ulb . ac . be/Nacr e/barre_database . htm 

3 ZAMS models are supposed to rotate like solid bodies, the cen- 
trifugal force is conservative and it is possible to define equipotential 
surfaces. This is no longer the case for shellular rotation. In that case 
however, the shape of the isobaric surfaces are given by the same ex- 
pression as the one giving the shape of the equipotentials for solid body 
rotation provided some changes of variables are performed (see Meynet 
& Maeder 1997). 
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Fig. 1. Variations of the radius as a function of initial mass, for vari- 
ous metallicities for non-rotating stars. The inset zoom shows the small 
difference of radii in the case of the 1 Mq models. 



where R e is the equatorial radius, one easily obtains from Eq. (1) 
that 



p \ fi 2 / 



P 



\ fi 2 ) 



(4) 



The above equation relates the inverse of the oblateness / to R p . 

As recalled above, the critical angular velocity, fi crit , corre- 
sponds to the angular velocity at the equator of the star such that 
the centrifugal force exactly balances the gravity. The classical 
critical angular velocity or the Q-limit (to distinguish it from the 
fir-limit as defined in Maeder & Meynet 2000a) in the frame of 
the Roche model is given by 



fieri t — 



2\' 2 



GM 

p,cnt 



(5) 



where /? p>cr j t is the polar radius when the surface rotates with the 
critical velocity. The critical velocity is given by 



2 GM V 



3R, 



p.cnt 



(6) 



Replacing fi in Eq. (4) by (fi/fi crrt ) • fi CI i t and using Eq. (5), one 
obtains a relation between R ? , f and R ?yCI n, 



Rr, — —R, 



p.crit 



Thus one has 



gcritf 3 (2{f-l)) l/3 



fi 

^crit 



3\ 3/2 //W\ 3/2 /2(/. 



R n 



1) 



P 



1/2 



(7) 



(8) 



With a good approximation (see below) one has that R p , C nt/Rp - 
1 and therefore 



fi 

ficrit 



3\ 3/2 /2(/ - 1)V /2 

~P~ 



(9) 



This equation is quite useful since it allows the determina- 
tion of fi cr it from quantities obtained with a model computed for 
any fi < fi crlt , however the closer to the critical value, the better. 
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It is not the case of Eq. (7) which involves /? PjCr it whose knowl- 
edge can only be obtained by computing models at the critical 
limit. Below we shall see to what extent the use of expression 
(9) provides a good approximation of O cr i t . 

Setting v the velocity at the equator, one has that 



1.05 



Q R„ R D R, 



p.cnt 



^crit ^c^e,crit ^crit ^p.crit ^e,crit 



(10) 



where ^ e , C rit i s tne equatorial radius when the surface rotates with 
the critical velocity. Using Eq. (7) above, and the fact that in the 
Roche model /? p , C rit/^e,crit = 2/3, one obtains 



— = 7?-2tf-l) 

Vcnl \ ^crit 



1/3 



(ID 



An interesting quantity is the ratio of the centrifugal acceler- 
ation a cen to the gravity g e at the equator 



^cen 



GM 



Q \ 2 (2\ 3 



R n 



R, 



p.crit 



(12) 



where we have used Eq. (5) and divided/multiplied by R 3 cril - 
Replacing /? p>c rit/^p by its expression deduced from Eq. (7), we 
obtain 



— =2(/-l). 



(13) 



We can check that at the critical limit, when / = 3/2, then a cen = 

ge- 

3. Basic parameters on the Zero Age Main 
Sequence 

We consider that a star has arrived on the ZAMS, when a mass 
fraction of 0.003 of hydrogen has been transformed into helium 
at the centre. At this stage, the star is supposed to have a solid 
body rotation. We present here the stellar radii, the total luminos- 
ity, the variation with the colatitude of the effective temperature 
and the moment of inertia of the models. We will then show the 
relations between the initial rotation and the angular momentum 
content, tWfcrit and Q ini /Q crit . 

3.1. Radii, luminosities and effective temperatures 

Considering first the case without rotation, we can see in Fig. 1 
that there is a relation between radius and mass of the form 
R ~ M° 6 at standard metallicity. This relation becomes shal- 
lower with decreasing Z, down to R ~ M° 3 at Z = 0. This comes 
mainly from two causes: 1) first, at lower metallicity, the opac- 
ities are lower; 2) the CNO elements are less abundant, or even 
absent in metal-free stars. In order to compensate for the energy 
radiated by the surface when smaller amounts of CNO elements 
are present, the central regions must extract energy from nuclear 
reactions at a higher temperature regime. To reach this regime, 
the star contracts more during the pre-MS phase and is thus more 
compact on the ZAMS. For example, in the case of the 60 Mq 
models, the radius at Z = is reduced by a factor 3.36 in com- 
parison with the radius at standard metallicity. The 1 Mq models 
are running only on pp-chains and show very small differences 
in their radii: only about 5% between Z = 0.020 and Z = (see 
the inset zoom in Fig. 1). 

In Fig. 2, we present the variation of the polar radius R p as a 
function of oj = Q/Q cr i t , normalized to its non-rotating value, for 
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Fig. 2. Variations of the polar radius as a function of the ratio w = 
il/n crit , normalized to the non-rotating value, for various masses at 
standard metallicity. 



the standard metallicity models. We see that increasing to tends 
to decrease the radius in the case of the 1 to 20 Mq, a trend which 
can be explained as follows: rotation plays the role of a sustain- 
ing force in the star. Thus, the central temperature and density 
in rotating stars are lower, decreasing subsequently the luminos- 
ity and the radius. The change of the polar radius for the 1 Mq 
model amounts to about 5%. In the case of the 3, 9 and 20 Mq 
models, the assumption of R p , cr i t /R p = 1 (Sect. 2.2) is absolutely 
valid, since the decrease amounts only to 1.5%. Why are nearly 
completely radiative stars as the 1 Mq stellar models more de- 
formed for a given value of co than stars with a significant con- 
vective core? It comes from the fact that the density gradient in 
convective zones are shallower than in radiative ones, allowing 
more mass to be contained in a given volume. This strengthens 
gravity and the resistance to deformation. 

For the 60 Mq models, the change of the polar radius is of 
the same amplitude as for the 1 Mq model but in the opposite 
direction, showing an increase of about 4.5%. How could be ex- 
plained such a difference of behaviour between the lowest and 
highest mass models? In the high mass star range, radiation pres- 
sure contributes more to the total pressure. It is the highest in the 
polar regions where the effective gravity (gravity decreased by 
the centrifugal acceleration) is the highest (von Zeipel theorem). 
A look at Fig. 6 shows that T e ff (pole) is higher when u is higher 
(10% between u = 0.1 and cj = 0.99), and therefore P md ~ T 4 
is stronger. This tends to inflate the polar radius in massive stars. 

Rotation deforms the star, the ratio R £ /R p increasing with o> 
as shown in Fig. 3. We have plotted the simplified relation ob- 
tained by Eq. (9), assuming R p , c dt/Rp = 1 (continuous line), as 
well as the relation obtained using Eq. (8) applied to two differ- 
ent masses (1 and 60 Mq) at Z = 0.02. In this case, we observe 
the direct effect of the behaviour seen in Fig. 2. For the 1 Mq, 
the decrease of R p with increasing oj makes / higher than the 
continuous line, joining it only at the critical limit since, in the 
Roche model frame, /? e ,crit/^p,crit is always equal to 1.5, while 
the inverse behaviour of the 60 Mq leads to the inverse result. 
Let us note also that in all cases, the equatorial radius becomes 
longer than 1.1 times the polar one only for o> > 0.75. 

In Fig. 4, we show the behaviour of the luminosity L(oS) of 
our 60 Mq models as a function of co for the different metallici- 
ties considered. The value of the luminosity is normalized to its 



Ekstrom et al.: Origin of Be stars 



5 



1.4 



^3 

a. 

Pi 



1.2 




<6 U 



0.2 



0.4 0.6 

"=o/o crit 



0.8 



Fig. 3. Variation of the ratio R e (cu)/R p (co) (equatorial over polar radius) 
as a function of <i> = Q/fl c . The continuous line shows the relation given 
by Eq. (9), assuming R v , C nt/R v = 1. The dashed lines show the relations 
for the Z = 0.020 modeis with 1 and 60 Mq using Eq. (8). 
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Fig. 4. Variations of the total luminosity as a function of oj, normalized 
to the non-rotating value, for 60 Mq at various metallicities. 
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Fig. 5. Variations of the total luminosity as a function of a>, normalized 
to the non-rotating value, for various masses at standard metallicity. 



4.6 



4.5 - 



E-. 
M 

O 4.4 



4.3 




20M r 



0.020 



30 



[°] 



60 



90 



Fig. 6. Variations of the effective temperature T e ff as a function of colat- 
itudinal angle 9, for the various values of rotational rates, in the 20 Mq 
model at standard metallicity. 



non-rotating value. One sees that the luminosity decreases by at 
most 5 to 10% when rotation increases up to the critical limit. As 
mentioned above, it is due to the fact that the centrifugal accel- 
eration helps in sustaining the star against the gravity and allows 
the luminosity to take a value corresponding to a non-rotating 
lower initial mass star. We can see that the effects of rotation are 
stronger at lower metallicity. This is related to the fact that low 
metallicity stars are more compact. The regions having a temper- 
ature high enough to undergo nuclear burning are more extended 
in the star and thus overlap regions that are more affected by ro- 
tation. 

Fig. 5 is the same as Fig. 4, but this time at fixed metallicity 
(Z = 0.020) and varying the mass. The 1 Mq is the most affected 
by rotation, with a decrease of 22.5% of its luminosity at break- 
up. The higher masses show a decrease of 4.5 to 7% only. This 
again is related to the fact that in the 1 Mq model, the energy 
production through pp-chains is more extended in the interior 
of the star. As above, a non-negligible part of the luminosity is 
produced in regions that may suffer from the effects of rotation. 

Figure 6 shows the variation with colatitude of the ef- 
fective temperature at the surface of a 20 Mq star for vari- 
ous rotation rates at standard metallicity. The ratio (r ei f(pole)- 



r ei f(equator))/r e ff (equator) becomes superior to 10% only for 
a> > 0.7. Near break-up, the effective temperature of the polar re- 
gion is about a factor two higher than that of the equatorial one. 
We may compare this result with the recent work of Espinosa 
Lara & Rieutord (2007) who have studied the structure and dy- 
namics of rapidly rotating stellar models with a two-dimension 
code. They have computed a model with a> = 0.82 for which the 
ratio (r e ff(pole)-r e ff(equator))/r e ff(equator) is 14%. Our closer 
model (co = 0.80) shows a slightly higher ratio of 16% though 
we should expect a smaller value due to the lower at. Let us note 
that the 2D simulation has been carried out within a spherical 
box which, according to the authors, may soften the latitudinal 
variations. It would be interesting to check how much of the dis- 
crepancy remains would the 2D computation be done without 
such a container. 

3.2. Inertia, angular momenta and rotational energies 

Fig. 7 presents the variations of the moment of inertia with the 
initial mass of the models. These variations show almost no de- 
pendence on o>, so only the case of o> — 0.1 is plotted. The mo- 
ments of inertia have been obtained by summing the contribu- 
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Fig. 7. Variations of the moment of inertia as a function of initial mass 
for various metallicities (models with Q/Q. ait = 0.1). 
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Fig. 8. Variations of the total angular momentum content of the models 
on the ZAMS as a function of initial rotation rate a> = Cl/Cl^t for var- 
ious metallicities, plotted according to the initial mass of the models. 
The coding for Z is the same as in Fig. 7 



tions of each shell ;' of the model: / = £ ( - (2/3) r 2 dm ; , where r, 
is the mean radius of the shell, i.e. the radius in P2( cos 8) = 0. 
Thus, r, accounts for the rotational distortion up to the second or- 
der Legendre polynomial, and the above way of calculating the 
moment of inertia is correct at this level of approximation. Let 
us mention that rotation deforms only the outermost layers of the 
model, which contain very little mass. This is why the moment 
of inertia does not change much with or without rotation. 

We note that, at a given metallicity, the moment of inertia / 
increases with the initial mass. The increase amounts to about 4 
orders of magnitudes between 1 and 60 Mq. As can be seen on 
Fig. 7, the relation between / and M is almost linear in logarith- 
mic scale, with a slope of 2.2 for the non-zero metallicities and 
of 1.7 for Z = 0. The difference of the slopes can be linked to 
the difference in the way the radii increase with increasing mass. 
From the mass-radius relations obtained in Sect. 3.1 (R ~ M - 6 
at standard metallicity and R ~ M° 3 at Z = 0) we indeed expect 
a relation / ~ M 2 2 and / ~ M l 6 at Z = 0.02 and respectively. 
For a given initial mass, the moments of inertia are smaller at 
lower metallicities, as expected due to smaller radii at lower Z. 

The variation of the total angular momentum content (X = 
2i h Q) as a function of oj is shown in Fig. 8. Note that the high- 
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Fig. 9. Variations of the ratio between rotational and potential energy 
as a function of the ratio a> = fJ/£2 crit for various metallicities, plotted 
according to the initial mass of the model. The coding for Z is the same 
as in Fig. 7 



est point on the right of the figures gives an estimate of the max- 
imum angular momentum content on the ZAMS (for solid body 
rotation). Since the moment of inertia of the star does not depend 
much on the rotation velocity, and since on the ZAMS, we sup- 
pose solid body rotation, one obtains quasi linear relations be- 
tween X. and u>. We note that the metallicity dependence remains 
modest at least in the range of metallicities between 0.00001 and 
0.020. In the 9 to 60 Mq range, the value of a> being kept con- 
stant, there is at first a slight increase of the angular momen- 
tum content when decreasing the metallicity (compare the curve 
with the black squares for the Z = 0.020 models with the curve 
with the triangles corresponding to 0.002 models). Then for still 
lower metallicities, the angular momentum decreases with the 
metallicity (the curve for the Z = 0.00001 models overlaps the 
Z = 0.020 models) reaching its lower values for the metal-free 
stellar models. 

In the mass range between 9 and 20 Mq, stars having the 
same angular momentum content on the ZAMS would have very 
similar value of co whatever their metallicity between 0.00001 
and 0.02. Only if the metallicity is zero, would the value of to 
corresponding to the same angular momentum content be much 
higher. As a numerical example, a 20 Mq on the ZAMS with £,- 
0.3T0 53 cm 2 g sec -1 has a ratio u> = 0.5 when Z is comprised 
between 0.00001 and 0.020. This angular momentum content 
corresponds tow= 0.8 when Z = 0. 

The rotational kinetic energy (E 10t = (l/2)/,Q 2 ) 
expressed as a fraction of the binding energy (^poti = 
—GYti (Mi/ri) d mi) is shown in Fig. 9. The rotational energy 
amounts to at most a percent of the binding energy. This is con- 
sistent with the well known fact that the effects of the centrifugal 
acceleration remains quite modest on the hydrostatic structure of 
the stellar interior, even when the surface rotates with a velocity 
near the critical one. For a 60 Mq stellar model at Z = 0.020 
the binding energy is 1.8 • 10 51 ergs, its Z = counterparts has 
a binding energy equal to 5.3 • 10 51 ergs, i.e. about a factor 3 
higher. Let us recall for comparison that the binding energy of 
a neutron star is of the order of 10 53 ergs, about two orders of 
magnitude greater. 
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Table 1. Angular velocity, rotational period, polar radius and equatorial 
velocity at the critical limit on the ZAMS (solid body rotation). 



Fig. 10. Relation between u/u crlt and to = fi/fi cr i t obtained in the 
frame of the Roche model. The continuous line is obtained assuming 
R p ,cnt/Rp = 1 (see text and Eqs. (9) and (11)). The dot-dashed lines 
show the relations for the Z = 0.020 models with 1 and 60 Mq using 
Eq. (8). The dotted line is the line of slope 1. 
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Fig. 11. Variations of the critical equatorial velocity on the ZAMS as a 
function of the initial mass, for various metallicities. 
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Fig. 12. Variations of the surface equatorial velocity on the ZAMS as a 
function of initial rotation rate w = Q./Cl ait for various metallicities and 
masses. The coding for Z is the same as in Fig. 11. 
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3.3. Relations between Q/Q cr i t and v/v C nt 

The relations between v/v cl i t and Q/fi C rit obtained in the frame 
of the Roche model (see Eq. (11)) for the 1 and 60 Mq stellar 
models at Z = 0.02 are shown in Fig. 10. In case we suppose that 
the polar radius remains constant (/? p>C rit/^p = 1), then Eq. (9) 
can be used and one obtains a unique relation between £2/f2 cr i t 
and f/fcrit, independent of the mass, metallicity and evolution- 
ary stage considered (continuous line). One sees that the values 
of u/fcrit are smaller than that of 0/Q C rit by at most ~25%. At 
the two extremes the ratios are of course equal. 



3.4. Velocities on the ZAMS: v crit , v 

Table 1 gives for the various models the angular velocity, the pe- 
riod, the polar radius and the equatorial velocity at the critical 
limit on the ZAMS. Let us note that for obtaining the value of 
Qcrit or of whatever quantities at the critical limit, it is necessary 
to compute models at the critical limit. Here we have computed 
models very near but not exactly at the critical limit, since at 
that point numerical singularities are encountered. As explained 
in Sect. 2.2, using Eq. (9) we could have estimated Q cr ; t from 
whatever model computed with a lower initial velocity. Doing 
this for the whole range of initial masses, metallicities and ve- 
locities, we would have obtained values within 10% of those 
shown in Table 1. For the models between 3 and 20 Mq, the er- 
ror would be always inferior to 3.5% because in this mass range, 
Rp(cu)/R p (0) remains very near 1 (see Fig. 2) and Eq. (9) is a 
good approximation of Eq. (8). For the 1 and 60 Mq the errors 
would be larger, amounting to a maximum of 9.3 and 6.7% re- 
spectively, showing the need of using rather Eq. (8) if co « 1. 
In general the error decreases drastically when o> increases, so 
in order to obtain the values given in Table 1, we have used our 
faster rotating model, deriving the value of £2 CI j t (using Eq. (9)) 
with an error of less than 0.001%. The critical period, P cr i t is then 
simply given by P cr i t = 2n/Q a i t . The polar radius is obtained 
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Fig. 13. Evolutionary tracks in the theoretical Hertzsprung-Russell diagrams for all the models computed beyond the ZAMS. The effective tem- 
perature corresponds to a surface averaged effective temperature. The computations are stopped when the star reaches the critical limit or before 
the beginning of the core He-burning phase. 



from Eq. (5) and the critical equatorial velocity is deduced from 

Went = 1-5 ^p,crit ^crit- 

The variation of the critical velocity as a function of the ini- 
tial mass and metallicity is shown in Fig. 1 1 . For a given initial 
mass, the critical velocity is higher at lower Z: e.g. for the 60 Mq 
model, the critical velocity goes from slightly more than 800 km 
s" 1 at Z = 0.020 to more than 1500 km s" 1 at Z = 0. It is a 
consequence of the smaller radii of stars at low metallicity (see 
Eq. (6)) as can be seen from Fig. 1. 

At a given metallicity, the critical velocity becomes higher 
for higher initial masses. This reflects the fact that the ratio of 
the mass to the radius increases with the mass of star (cf. Sect. 
3.1). The critical rotation periods are between 2.4 and 25 hours 
for the whole domain of masses and metallicities explored in the 
present work. Let us just mention here that, as will be described 
in Sect. 5, the critical velocity evolves with time, so the values 
given here are values at the beginning of the star's evolution and 
should not be used throughout the whole MS phase. 

Fig. 12 shows relations between v and a> on the ZAMS for 
different initial mass models at various metallicities. The cor- 
responding values are also given in Table 2. We see that, for a 
given value of at, higher surface velocities are obtained at lower 
metallicities and for higher initial masses. 



4. Evolution on the Main Sequence 

4.1. Hertzsprung-Russell diagram and lifetimes 

Evolutionary tracks in the Hertzsprung-Russell diagram for the 
present models are shown in Fig. 13. The models have been 
evolved until either the reaching of O cr ; t or just before the onset 
of central Helium burning. For each mass, the effect of reducing 
the metal content shifts the tracks towards the blue part of the di- 
agram and, for 3 and 9 Mq, towards higher luminosities. In each 
metallicity subgroup, the effect of raising up the rotational rate 
shifts the tracks towards lower luminosities and cooler temper- 
atures on the ZAMS. When the star evolves, rotational mixing 
continuously refuels the convective core in fresh hydrogen and 
brings newly produced helium in the radiative zone. The enrich- 
ment of the radiative envelope in helium makes the star brighter 
and bluer by decreasing the opacity. If we look at what happens 
to the convective core (Fig. 14), we see that at first, a higher ro- 
tational rate leads to a smaller convective core, due to the lower 
central temperature. Then, we see that the decrease of the core is 
slower because of the refuelling in hydrogen. At the end of the 
MS, the core is much larger in the fast rotating models. In that 
respect rotation acts as a core overshoot. 

In Table 3 the lifetimes of our stellar models are given as 
well as the time at which the critical limit (if reached) is first 
encountered during the MS phase. For the models which were 
not pursued until the end of the MS phase, we give an estimate 
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Table 2. The first line indicates for each initial mass and metallicity, the 
equatorial velocity in km s -1 on the ZAMS. For the 3, 20 and 60 Mq 
stellar mass models the second line gives the time-averaged value of the 
velocity over the adjustment period from solid to non-solid rotation at 
the beginning of the MS phase (see Sect. 4.2). 
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of the MS lifetime (this is indicated by an asterisk) obtained by 
inspections of the variation of the central hydrogen mass fraction 
with time. Note that we suppose that the star, once it reaches the 
critical limit, remains in the vicinity of this limit. This behaviour 
has been confirmed by previous computations in which models 
were pursued beyond the point where they first reach the critical 
limit (see Ekstrom et al. 2006). Most of the 60 Mq stellar models 
encounter the critical limit very early. In that case only the time 
of reaching the limit is given. We see that for u> = 0.7 the MS 
lifetime is increased by 15 - 23% (depending on the initial mass 
and metallicity) in comparison with the models with u = 0.1. 

4.2. Evolution of the equatorial velocity 

The evolution of the rotation velocities at the stellar sur- 
face depends of course directly on the movements of expan- 
sion/contraction of the star. But besides this, it also strongly de- 
pends on two factors: the internal mechanism of angular momen- 
tum transport (internal coupling) and the mass loss. 

An extreme case of internal angular momentum transport is 
the one which imposes solid body rotation at each time in the 



0.5 




age [Myr] 



Fig. 14. Evolution of the size of the convective core during MS, for the 
various values of rotational rate, in the 20 Mq at standard metallicity. 

course of the evolution of the star. A strong coupling is then 
realized between the contracting core and the expanding enve- 
lope. The spinning up of the central regions is transmitted to the 
stellar surface. Thus, during the MS phase, the surface velocity 
approaches the critical limit and may even reach it, depending 
on the initial conditions (see Sackmann & Anand 1970; Langer 
1997; Meynet et al. 2007). Another extreme case is the case of 
no transport of angular momentum. Each stellar layer keeps its 
own angular momentum. The variation of Q is then simply gov- 
erned by the local conservation of the angular momentum and, at 
the surface, Q. ~ R 2 . When the radius increases, the surface an- 
gular velocity decreases more rapidly than the classical critical 
angular velocity (see Eq. (5)), thus the star evolves away from 
the critical limit. In the present models, without magnetic fields, 
the situation is intermediate between these two extreme cases. A 
moderate coupling is exerted mainly by meridional circulation, 
which is more efficient than shear transport, as far as transport 
of angular momentum is concerned (see e.g. Meynet & Maeder 
2000). 

The second important effect governing the evolution of an- 
gular momentum is mass loss. The strong mass losses during the 
MS phase, expected for high mass stars at high metallicities, de- 
crease the angular momentum of the stars. The importance of 
the decrease of the angular momentum depends on how matter 
is ejected. In the extreme case where all the matter would be 
ejected along narrow jets centred on the rotational axis, nearly 
no angular momentum would be lost, while if all the matter is 
ejected near the equator plane, this maximizes the amount of an- 
gular momentum lost. Depending on the geometry of the wind, 
the surface velocity may thus decrease or increase. In the present 
work, we did not account for the effects of wind anisotropies 
(which would be important only for a relatively narrow range of 
masses and velocities, see Meynet & Maeder 2003), thus, in the 
present case, strong mass loss always leads to a decrease of the 
surface velocity. 

Figure 15 shows the evolution of the surface velocity ex- 
pressed as a fraction of the critical velocity during the MS phase 
for 20 Mq stellar models. Quite generally the curves can be de- 
composed into three parts, whose importance depends on the ini- 
tial rotation and metallicity: 

1. At the very beginning, there is a short adjustment period, 
which lasts for a few percent of the MS lifetime, during 
which meridional circulation transports angular momentum 
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Table 3. Lifetimes (in Myr) during the MS phase. An asterisk means 
that the value was estimated (see text). 



z n/n c 





0. 1 


0.3 


0.5 


0.7 


0.8 


0.9 


0.99 










3M© 











215.6 


227.8 


237.8 


249* 


255* 


264* 


275* 










(223.7) 


(206.2) 


(189.0) 


(170.6) 


1(T 5 


224.9 


243.6 




261.5 


265.6 


271* 


271* 














(270.6) 


(213.4) 


0.002 


274.2 


301.3 


312.3 


320.9 


325.4 


331* 


331* 














(316.4) 


(237.9) 


0.020 


357.0 


399.6 


424.4 


439.2 


444.4 


454* 


484* 














(446 7) 


(302 9) 





















20.31 


22.11 


23.24 


24.31 


24.74 


25* 


25* 














(24.77) 


(16.12) 


to- 5 


26.62 


28.70 


29.98 


30.85 


31.29 


31.8* 


32.4* 














(31.30) 


(25.03) 


0.002 


27.77 


30.11 


31.32 


32.26 


32.87 


33.5* 


34.0* 














(31.84) 


(26.00) 


0.020 


27.59 


30.33 


31.62 


32.62 


33.0* 


34.0* 


34.0* 












f39 Ql \ 


(28 89) 


(23 05) 










20M Q 











8.003 


8.554 


8.868 


9.199 


9.324 


9.434 


9.44* 
















(8.616) 


to- 5 


9.099 


9.853 


10.213 


10.54* 


10.68* 


10.73* 


10.73* 










(10.533) 


(9.626) 


(7.843) 


(5.651) 


0.002 


9.303 


10.117 


10.527 


10.85* 


11.00* 


11.06* 


11.15* 










(10.840) 


(10.033) 


(8.444) 


(6.712) 


0.020 


8.663 


9.521 


9.928 


10.39* 


10.66* 


10.96* 


10.96* 










(10.389) 


(10.662) 


(10.937) 


(7.892) 










6OM 











3.633 


3.802 


3.924 


















(3.386) 


(2.538) 


(2.259) 


(1.173) 


to- 5 


3.874 


4.112 


















(3.696) 


(2.746) 


(2.200) 


(1.525) 


(0.783) 


0.002 


3.994 


4.251 


















(4.432) 


(2.818) 


(2.336) 


(1.717) 


(1.087) 


0.020 


3.714 


3.978 


















(3.960) 


(2.426) 


(2.031) 


(1.595) 


(1.179) 



from the outer parts of the star to the inner ones. This slows 
down the surface of the star. Then, in the interior, shear tur- 
bulence becomes active and erodes the gradients built by 
the meridional circulation. Under the influence of these two 
counteracting effects the Q-profile converges toward an equi- 
librium configuration (Zahn 1992; Denissenkov et al. 1999; 
Meynet & Maeder 2000). From Fig. 15, one sees that for 
lower initial values of v/w cr ; t the decrease is much weaker 
than for higher values. It may even completely disappear. It 
illustrates the dependence of the meridional velocity on O. 
This feature also shows that, in the present theoretical con- 
text, we cannot have stars which would remain near the crit- 
ical limit starting from the ZAMS. Indeed, our most extreme 
case, the model with £l/Q CI \t = 0.99 evolves away from the 
critical limit just after the ZAMS. Of course we can won- 
der to what extent this behaviour is dependent on the chosen 
initial conditions. For instance would this behaviour be the 
same, if one would have computed realistic pre-MS evolu- 
tion? Probably in that case the model would not have solid 
body rotation on the ZAMS and thus the above feature is 
in part artificial. It is however difficult to say more due to 
the complexity of the processes affecting the angular mo- 




Time [Myr] Time [Myr] 

Fig. 15. Evolution of the ratio v/v a n during the core hydrogen-burning 
phase for the 20 Mq models at various metallicities. The computations 
were stopped as soon as the star reaches the critical limit or just before 
the beginning of the core He-burning phase. 



mentum of stars during their formation process. The point 
which has to be kept in mind is that in the present models 
we cannot have non-evolved critically rotating stars. Let us 
stress that these stars should be anyway very rare and their 
absence should not have any impact on the frequency of stars 
predicted to rotate at or near the critical velocity. The above 
effect implies also that the velocities on the ZAMS given in 
Table 2 are probably quantities not directly comparable to 
observed ones. Indeed, in the most optimistic case, where 
ZAMS models would indeed rotate as solid body, the chance 
to see stars at the very beginning of their relaxation period 
is very small. More realistic values of the initial velocities 
would correspond to time averaged values computed over the 
period during which about ten percent of the initial content 
in hydrogen has been burnt in the core. These values are also 
given in Table 2. It is important to note that models treating 
advection as a diffusion would miss this first relaxation pe- 
riod. Indeed during this phase, meridional currents build the 
gradient of Q., a process that diffusion cannot do, diffusion 
always eroding any gradient. 
2. As explained for instance in Maeder & Meynet (2000b), af- 
ter the first adjustment period, a large outer cell of meridional 
circulation sets in, transporting angular momentum from the 
inner parts of the star to the outer ones 4 . If this transport 
is rapid enough and the stellar winds not too intense, then 
v/Vcnt increases (an important distinction has to be made 
here: it is the ratio v/v CI n which increases, not necessarily 
the equatorial velocity itself, see Fig. 18 and Sect. 5 for dis- 
cussion on this point). It is the case for all the models with 
to > 0.3 at Z > 0.00001. For metal-free stars, the increase 
of y/v cr it only occurs for the models with a> > 0.7. This is 
a consequence of the Gratton-Opik effect (see the term in 
brackets in Eq. (14) making the meridional circulation ve- 
locity smaller at lower Z. 



4 In our model, during the MS phase, we have two cells of merid- 
ional circulation, an inner one where matter moves clockwise, and an 
outer one where matter moves counter-clockwise. The velocities in the 
outer cell are much larger than the velocities of the inner cell. As a con- 
sequence, the outer cell clearly dominates the transport of the angular 
momentum in the star, transporting it from inner regions towards outer 
ones. 
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Fig. 16. Comparison of the evolution of v/v aii (top) and of log M (bot- 
tom) for the 20 Mq models (left) and the 60 Mq models (right) at stan- 
dard metallicity, when using Kudritzki & Puis 2000 (continuous lines) 
or Vink et al. 2000 (dotted lines). 



3. At the end of the MS phase, the increase of vjv cn t acceler- 
ates when the star contracts after the blue hook in the HR 
diagram. The rapid change of the velocity at that stage is due 
to the rapid change of the radius, not to any transport mech- 
anisms. 

Looking at the right bottom panel of Fig. 15 (models at 
Z = 0.020), another interesting effect can be discussed. We see 
for instance that the 20 Mq stellar models with initial value of 
v I Vent between 0.5 and 0.7 (cj between 0.7 and 0.9) show a bump 
around an age of 9 Myr. This occurs because at that time, for 
these models, mass loss increases enough to remove efficiently 
the angular momentum brought to the surface by the meridional 
currents. The surface velocity thus evolves away from the critical 
limit, reaching it again only at the very end of the MS phase. This 
feature is not seen at lower metallicity because mass loss rates 
are lower. This trend was already put in evidence by Maeder & 
Meynet (2001). For a given degree of coupling, the mass loss 
rates play a most critical role in the evolution of the surface ro- 
tation. 

In order to check the validity of our results, we have com- 
puted another set of 20 and 60 Mq models at Z = 0.020 with the 
mass loss prescription of Vink et al. 2000 (hereafter Vink et al. 
2000). In the case of the 20 Mq models, the mass loss rates ob- 
tained in the new set are slightly lower than the ones obtained 
with Kudritzki & Puis 2000 (KP00). The a> > 0.8 models meet 
thus the break-up limit earlier in their evolution. The a> = 0.7 
model shows also a bump in the velocity, as described above, 
but then it stays far from break-up limit, even at the end of the 
MS unlike what was the case with the KP00 rates. This has for 
consequence that the minimum a> required to reach the break-up 
limit during evolution will be slightly higher with VdKLOO rates, 
while the fast rotators will spend a little more time at break-up 
limit (Fig. 16, left panels). But the qualitative result is generally 
similar to the one obtained with KP00 rates. In the 60 Mq case, 
the opposite happens: VdKLOO's rates are much stronger (by a 
factor 2 during the MS), and thus no model reaches the break-up 
limit anymore, even the one starting its evolution with w = 0.99 
(Fig. 16, right panels). The most massive stars' evolution is dra- 
matically influenced by the choice of the mass loss rates and it 
can be concluded that in that case, the comparison between the- 



Table 4. Time averaged equatorial velocities in km s -1 corresponding to 
various ZAMS values of oj = Q/Q. a n. The average is taken either on the 
whole MS phase or on the part of the evolution before the critical limit 
is reached. For these latter models, the mass fraction of hydrogen at the 
centre when the surface velocity reaches the critical limit is indicated in 
parenthesis. 



M 


Z 


0.1 


0.3 


w 
0.5 


= 

0.7 


n/n cril 

0.8 


0.9 


0.99 


3 





38 


108 


198 


294 


333 


373 


423 












(5e-2) 


(0.1) 


(0.17) 


(0.26) 


3 


10- 5 


33 


95 


167 


247 


297 


353 


416 
















(7e-3) 


(0.30) 


3 


0.002 


27 


80 


139 


207 


246 


298 


352 
















(0.1) 


(0.38) 


3 


0.020 


23 


69 


120 


181 


215 


258 


311 
















(3e-2) 


(0.38) 


9 





56 


151 


265 


394 


477 


569 


666 
















(le-5) 


(0.35) 


9 


10- 5 


39 


117 


204 


305 


363 


435 


515 
















(3e-2) 


(0.29) 


9 


0.002 


33 


99 


173 


256 


305 


369 


425 














(3e-4) 


(9e-2) 


(0.31) 


9 


0.020 


30 


90 


157 


236 


284 


341 


387 














(2e-2) 


(0.21) 


(0.36) 


20 





54 


169 


295 


446 


507 


588 


721 


















(0.13) 


20 


10- 5 


45 


144 


254 


389 


469 


557 


651 












(3e-3) 


(0.15) 


(0.33) 


(0.48) 


20 


0.002 


37 


120 


213 


323 


387 


456 


528 












(2e-3) 


(0.14) 


(0.30) 


(0.43) 


20 


0.020 


32 


107 


189 


279 


325 


369 


447 












(2e-4) 


(le-4) 


(6e-4) 


(0.30) 


60 





63 


248 


456 


645 


830 


930 


1167 












(0.20) 


(0.38) 


(0.45) 


(0.61) 


60 


10- 5 


63 


221 


408 


587 


687 


797 


890 










(0.16) 


(0.37) 


(0.47) 


(0.57) 


(0.67) 


60 


0.002 


48 


179 


323 


474 


552 


641 


727 










(5e-5) 


(0.40) 


(0.45) 


(0.54) 


(0.63) 


60 


0.020 


37 


140 


258 


400 


472 


547 


609 










(8e-2) 


(0.38) 


(0.44) 


(0.51) 


(0.56) 



oretical models and observations constitutes rather a check for 
the mass loss rates than anything else. 

Coming back to our initial set of models, we present in 
Table 4 the time averaged velocities during the MS phase or dur- 
ing the period before the reaching of the critical limit. This ta- 
ble is useful in order to connect the values of oj on the ZAMS 
with observable rotational properties during the MS phase (see 
Sect. 5). 



4.3. Surface enrichment 

In Fig. 17 the tracks are plotted in the plane (N/C)/(N/C) ini 
versus P, where P is the rotational period in hours. During the 
evolution the surface is progressively enriched in CNO burning 
products, i.e. enriched in nitrogen and depleted in carbon. At the 
same time, the rotational period increases. The hatched regions 
corresponds to zones where stars would rotate at or faster than 
the critical limit. We can note the following features: 

- As already mentioned in previous works (see e.g. Maeder & 
Meynet 2001), for a given value of the initial mass, of u> and 
for a given evolutionary stage, the surface enrichments are 
stronger in metal poor stars. This is related mainly to the fact 
that the gradients of Q are steeper in metal poor stars due the 
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Fig. 17. Evolutionary tracks in the plane surface N/C ratio, normalized to its initial value, versus the rotational period in hours for different initial 
mass stars, various initial velocities and for the metallicities Z=0.02 and 0.002. Positions of some periods in days are indicated at the bottom of 
the figure. The dotted tracks never reach the critical limit during the MS phase. The short dashed tracks reach the critical limit during the MS 
phase. The dividing line between the shaded and non-shaded areas corresponds to the entrance into the stage when the star is at the critical limit 
during the MS evolution. Big circles along some tracks indicate the stage when v/v ctit becomes superior to 0.7. If Be stars are rotating at velocities 
superior to 70% of the critical velocity, present models would predict that they would lie in the region comprised between the big circles and the 
dividing line. Beware the different vertical scales used when comparing similar masses at different metallicities. 



Ekstrom et al.: Origin of Be stars 



13 



600 



m 500 

6 



? °400 



300 



200l=-L 



v/v, 

y 



crit J 




0.8 



0.6 



0.4 



0.2 



4 6 
Time [Myr] 



Fig. 18. Evolution of the equatorial velocity v (continuous line), the crit- 
ical velocity u crit (dotted line) and their ratio (dashed line) for the 20 Mq 
at standard metallicity with to = 0.70 initially (u/u cr i t = 0.50). 



slower meridional currents in more compact stars. We note 
that surface enrichments for 3 Mq stellar models at Z = 0.02 
remain in the range of a few percents, even for high rotation 
rates. This means that for these low mass stars, observations 
will probably not be able to see such small enrichments. The 
same initial mass models at Z = 0.002 present enhancements 
of up to 45%. 

- Also the surface enrichments are higher in more massive 
stars, a feature which is due to higher efficiency of the shear 
mixing when the initial mass increases. 

- Stars with higher initial velocities reach higher N/C ratios. 
Let us recall here that the computations were stopped when 
the critical limit is reached, this explains why for instance, 
the maximum surface enrichment indicated for the 9 Mq 
stellar model at Z = 0.02 is lower for Q/Q cr it = 0.99 than for 
n/Q crit = 0.8 and 0.9. Would the model with Q/Q crit = 0.99 
be continued, it would reach higher values than those shown 
in Fig. 17. 

- In general, stars near the critical limit show nitrogen en- 
hancements. There are however cases where no or very small 
enhancements would be expected. These cases are: the stars 
starting their evolution with very high rotation rates and still 
being nearly non-evolved. For instance, in case of the 9 M 
stellar model at Z = 0.02 with Q/n ini = 0.99, v/v crit be- 
comes superior to 0.7 while still no significant surface en- 
richment in nitrogen has occurred. Another case where no 
observable surface enrichment is expected is for stars at 
Z = 0.02 and with initial masses around 3 M (see lower 
left panel). 

- No stars with masses between 3 and 20 Mq and with rotation 
periods longer than about 5 days are expected to be at the 
critical limit. 



5. Initial conditions for stars reaching near critical 
velocities 

Let us first recall that an important ingredient of the ratio v/Vcnt 
is the behaviour of u cr ; t itself. Fig. 18 shows how it evolves in 
time in the case of the 20 Mq model at standard metallicity 
with o) = 0.7 initially (v/v ct i t = 0.5). During evolution, u cr ; t 
gets lower, and in the present case the final value amounts to 
only half the initial one. It can be explained by the modifica- 
tions of the stellar parameters: some mass is lost and the radius 
is slowly inflating, so the surface gravity gets lower and can be 
more easily counterbalanced by the centrifugal force. In the ex- 
ample shown in Fig. 18, we see that even though the equato- 
rial velocity is slowly decreasing during the MS, the ratio v/v cr i t 
increases steadily. It is important to take into account this evo- 
lutionary effect when one wants to determine the v/v CT i t ratio 
at a given time. Keeping the initial value would lead to a com- 
pletely wrong result. Let us mention that in that perspective, our 
results are not contradictory with the observations of Huang & 
Gies (2006b) or Abt (2003) who find a clear spin down with 
evolution for the stars in their sample. Except the most extreme 
ones with a> = 0.99 - 0.90, all our standard metallicity models in 
the mass range 3 - 9Mq show a steady or decreasing equatorial 
velocity, though their f/f cr it ratios increase. 

From our computations, we can determine for each model 
specified by its initial mass and metallicity, the values of £2/f2 cr i t 
on the ZAMS required to reach given limits of v/v cr i t during the 
MS phase. These values are plotted in Fig. 19. Using Table 2, 
it is possible to translate these limits on Q/f2 cr i t into limits on 
the equatorial velocities on the ZAMS or at the beginning of the 
MS phase. Using Table 4 a similar correspondence can be ob- 
tained with values of the averaged velocity during the MS phase 
or during the phase before the reaching of the critical velocity. 

At the standard metallicity, we note that in general the values 
of Q/O cr it needed on the ZAMS to reach given values of f/i/ cr it 
decrease when the mass increases. This comes from the fact that 
meridional currents are more rapid in the less dense envelope 
of more massive stars. Let us recall that in the outer layer, the 
dominant term in the expression for the radial component of the 
meridional circulation U (r) is mainly (Gratton-Opik term) 



U(r) 



1 



Q 2 



2nGp 



n 



2,3 



GM 



(14) 



where p means the average over the considered equipotential. 
Thus the outer cell of the meridional circulation will be more 
active in less dense envelopes. In Fig. 19 upper left panel, the 
bump around the 20 Mq stellar model comes from the fact that 
this model suffers significant mass loss at the end of its MS evo- 
lution. This disfavours the reaching of near critical values (ex- 
cept at the very end of the MS evolution, see Fig. 15). In the 
case of the 60 Mq model, when the mass loss rates of KP00 
are used, the near critical values occur well before the star un- 
dergo similar enhancements of the mass loss rate and thus, the 
minimum values of f2/Q cr ; t are lower than for the 20 Mq. The 
situation is quite different if the mass loss rates from VdKLOO 
are used (see Fig. 16). 

For metallicities Z = 0.002 and 0.00001, we find that the 
minimum values of 0/Q cr i t required to reach near critical veloc- 
ities also decrease when the initial mass increases. More massive 
stars thus reach more easily the critical limit. We note also that 
for the 60 Mq at Z = 0.002, the minimum value of Q/fi cr it to 
reach the critical velocity is around 0.7, while at Z = 0.00001 it 
is lowered to values equal to 0.5. This is mainly due to the weak 
radiative winds of the metal-poor models. Since radiative winds 
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Fig. 19. Minimum values of u> = fi/£2 crit on the ZAMS required for reaching during the MS phase values of v/v al i superior to 0.4 (bottom line), 
0.5, 0.6, 0.7 (dashed), 0.8, 0.9 and 0.99 (top heavy line) respectively. The initial metallicity considered is indicated on each panel. 



are triggered by metal lines, a low metal content implies low 
mass loss rates. We have already seen in Sect. 4.2 that there is 
a competition between the efficiency of the mass loss (removing 
some angular momentum and thus decelerating the surface) and 
the efficiency of the meridional circulation (bringing the angular 
momentum from the core to the surface and thus accelerating the 
surface). From our results, we can say that even though the effi- 
ciency of the meridional circulation is lower at low Z, the effect 
of lowering the mass loss rates plays the major role and the total 
budget is in favour of a significant acceleration of the surface, at 
least down to Z = 0.00001. Let us note that at some points, the 
curves are superposed: for instance for the 9 Mq atZ = 0.002, 
the curve for reaching the critical limit, and those for reaching 
values of v/v cr \ t > 0.8 and 0.9 go through the same point. This 
simply reflects the fact that when such a star reaches the values 
f /fcrit ^ 0.8, it then reaches also the critical limit. 

For metal-free stars, one notes that for M > 9 M Q , the mini- 
mum values of Q/£2 cr j t to approach critical velocities during the 
MS evolution are in general higher than at higher metallicities. 
This reflects the greater compactness of these stars. The metal- 
free 3 Mq models on the contrary reach much more easily the 
critical limit than the models at higher metallicities. We can see 
that their tracks in the HR diagram (see Fig. 13) evolve quite 
blueward. This evolution is due to the fact that metal-free stars 



in this mass domain must contract over a great part of their MS 
evolution in order to compensate for the lack of CNO elements. 
In such a situation, the evolution of Q is mainly governed by 
the local conservation of the angular momentum (nearly no cou- 
pling) which favours the approach of the critical limit when the 
star contracts. Typically for a 3 Mq stellar model, the minimum 
value of Q/O cr it is around 0.6 at Z = for reaching the critical 
limit during the MS phase, while it is higher than 0.8 at standard 
metallicity. 



6. Rotation of stars in clusters of various ages and 
metallicities 

Let us call n(t, w) the number of MS stars in a cluster of age 
t having surface rotational velocities such that v/v cr n > w. We 
only consider stars brighter than about two magnitudes below the 
turn off and call N(t) the total number of MS stars in the same 
cluster, in the same interval of magnitudes. The ratio n(t, w)/N(t) 
can be expressed in the following way 

rM„ m , _ 

n(t,w) _ } Umin n(M)p(M,t,w)dM 
N(t) £~"ri(M)f{M,fm 
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Fig. 20. Left panel: distribution of the equatorial velocities as given by Huang & Gies (2006a) (continuous line). The dotted line shows the case 
when it is assumed that above 300 km s -1 , the velocities are underestimated (see text). The dashed line shows a flat distribution. At the bottom 
of the figure, the correspondence between ZAMS values of O / f2 crit and average velocities during the MS phase as deduced from our 9 Mq at 
standard metallicity (Table 4) is indicated. Right panel: cumulative distribution function of w = Q./Q. ait deduced from the data by Huang & Gies 
(2006a) with the same assertions as the left panel. 



where M m ,„ is the mass of stars which are at two magnitudes 
below the turn off, M max the mass of stars which are at the turn 
off, n(M), the number of stars born with an initial mass M, and 
p(M, t, w) the fraction of these stars which, at the age f, are on 
the MS and have surface velocities such that v / v all > w, f(M, f) 
is the fraction of stars of initial mass M which are on the MS at 
the time t. In the following we explain how these different values 
are obtained from our rotating stellar models. 



- M max : In the case of non-rotating stellar models, at a given 
metallicity we have a one-to-one relation between the initial 
mass and the MS lifetime. It suffices then to take the ini- 
tial mass value of the star which has a MS lifetime equal 
to the considered age of the cluster in order to obtain the 
mass of stars at the turn off, i.e. M max . When rotation is ac- 
counted for, the situation is more complicated, since the MS 
lifetimes increase with the initial velocities. Thus the maxi- 
mum mass of a star which can still be on the MS at a given 
age corresponds to the mass of the rotating star, having the 
maximum initial rotation velocity (here Cl/£l CI n = 0.99 on 
the ZAMS) and whose MS lifetime corresponds to the age 
considered. Of course in a given cluster there is little chance 
to have at turn off such a star due to the very small number 
of stars with such a high velocity. This effect will be taken 
into account through the initial velocity distribution function 
in f(M, t) (see below). 

- M m { n : To estimate the lower bound M m ; n , we have taken 
about 70% of the value of M max . We have indeed checked us- 
ing theoretical isochrones (Meynet et al. 1993) that for ages 
between 8 and 316 Myr, the typical mass of stars two visual 
magnitudes below the turn off is in that range. 



- n(M): The number of stars born in the mass interval between 

M min and M max is proportional to f^ ,mx M~ 235 dM, when a 
Salpeter's slope of the initial mass function is considered. 

- f(M, t): as explained above, stars with a mass M max will still 
be in their MS phase, provided they started their evolution 
with 0/n cr i t > 0.99. Stars of this mass but with slower initial 
rotation will have already evolved away from the MS. Thus 
the fraction of stars of mass M max which are still on the MS 
at time t is equal to 

^6J=1.00 

/W max , f) = I d p (uj)du) 

Jw min (M,0=0.99 

where <w m j n (M, f) is the minimum value of a> for stars of ini- 
tial mass M to be still on the MS at the time f, and d p the 
probability density for a star to start its evolution with a given 
value lo (See Fig. 20). For each mass and age, one can de- 
termine a> m i n (M, f), and similar integrals as the one shown 
above can be performed. 

- p(M, t, w): From figures similar to Fig. 15, we can extract for 
each initial mass, metallicity and rotation, the time at which 
w/vcrit becomes superior to w (note that once v/v cr \ t > w, it 
remains superior to it for the rest of the MS phase 5 , the only 
exception being the 20 Mq at Z = 0.020). From such data, 
one can deduce that for a star of 9 Mq with Z = 0.020 and an 
age equal to 25 Myr, the minimum value of f2/Q cr ; t required 
on the ZAMS to reach values of v/v CI \ t superior to 0.7 during 
the MS phase is 0.83. Using an initial distribution of Q/£l cr it 
(see Fig. 20), one can deduce the value of p(M, t, w). 



5 Note that we used the extrapolated lifetimes (marked with an aster- 
isk in Table 3), for the models whose computation was stopped before 
the completion of the MS phase. 
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Fig. 21. Variation of the fraction of stars with vjv aA superior or equal to 0.4, 0.5, 0.6 (three thin continuous lines from the top), 0.7 (thick dashed 
line), 0.8, 0.9 (two thin continuous lines) and at the critical limit (bottom thick solid line) in clusters of various ages and metallicities. At the bottom 
of each panel, the range of MS lifetimes for the 3, 9 and 20 Mq stellar models are indicated. 



We need now to specify the function we use for the initial 
distribution of £2/f2 cl i t . Huang & Gies (2006a) provide from their 
observations of 496 OB-type stars a distribution of equatorial ve- 
locities corrected for the sin i effect. This distribution is shown 
in Fig. 20 (left). Of course this distribution does not correspond 
to a distribution on the ZAMS, stars being picked up at different 
phases of their evolution during the MS phase. We can however 
reasonably consider that these equatorial velocities are represen- 
tative of the values of the average velocities of stars in the MS 
phase. We can then associate a value of a> on the ZAMS to each 
these values using the relation obtained from our 9 Mq at stan- 
dard metallicity between the averaged equatorial velocity during 
the MS phase and a> on the ZAMS. Proper normalization have 
been applied to transform the probability densities for velocities 
into probability densities for a>, i.e. the normalization has been 

D p (v)dv = J () d p (cj)dco = 1, where D p and 
d p are the probability densities for respectively v (see Fig. 20, 
left panel) and co. The cumulative distribution of a> on the ZAMS 
we obtained in this way is shown in the right panel of Fig. 20. 

In Fig. 21 the variation of the number fraction of stars at the 
critical limit (continuous thick lines) in clusters of different ages 
is shown. We have plotted also the fraction of stars having v/v CI i t 
values superior to given limits (see caption). Due to the coarse 
grid of mass we used here, we shall not attach too much impor- 
tance to the details of these curves but will comment on their 
general behaviour. A few interesting features can be underlined: 



- At standard metallicity, we note a bump around ages 20 - 25 
Myr for the curves corresponding to the critical limit and 
for f/f C rit ^ 0.9 and 0.8. This kind of behaviour can be 
explained in the following way: at younger ages than the 
bump, mass loss disfavours the reaching of the critical limit, 
at older ages than the bump, meridional currents are not effi- 
cient enough for transporting outwards the angular momen- 
tum (more dense envelopes in stars of lower initial masses). 

- For metallicities Z = 0.002 and 0.00001, we note that the 
bump shifts to younger ages and reaches higher values. This 
is due to the smaller mass loss rates experienced by stars at 
low metallicity. This allows more massive stars, which have 
less dense envelopes and thus stronger meridional currents, 
to undergo a more rapid acceleration of their surface and to 
reach more easily the critical limit. 

- We note that the curves are very similar at the metallicities 
Z = 0.002 and 0.00001. Thus, from the point of view of 
the surface velocities, we would not expect any significant 
changes over this metallicity range. 

- For metal-free stars the curves are very different than at 
Z = 0.00001 and 0.020: at ages younger than about 20 - 32 
Myr, stars have more difficulty to reach the critical limit. This 
comes from the fact that the absence of metals make these 
stars very compact, with dense outer layers and very slow 
meridional circulation currents (note that would the mass 
range be extended towards higher masses, it is possible that 
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this bump would appear at very young ages). At ages older 
than about 80 Myr, on the other hand, stars at the critical 
limit are much more numerous than at higher metallicities. 
This comes from the fact that stars with initial masses equal 
or below 3-5 Mq contract for a great part of their MS phases 
(see the metal-free 3 Mq stellar models) and thus can reach 
quite early the break-up limit. 
- At Z = 0.020, the fraction of stars at the critical limit is 
of the order of at most 5% and can reach values of slightly 
more than 10% at lower metallicities. The number fraction 
of stars with v/v cr i t > 0.9 tightly follows the curve of stars at 
the critical limit, indicating that when this value is reached, 
the surface velocity rapidly evolves towards the critical limit. 
The fraction of stars with v/v ct i t > 0.7 reaches values up to 
20-25% and this at an earlier age: around 10 Myr and below. 

7. Discussion of the results 

7.1. Sensitivity on the ZAMS model 

We recalled in Sect. 4.2 that our way of building our ZAMS 
models makes it impossible to maintain a high value of v/v CI \ t 
from the ZAMS all along the beginning of the MS phase. The ac- 
tivation of the meridional circulation imposes a rapid slow down 
of the surface and thus, the star needs to somewhat evolve in 
the MS phase in order to approach again the critical limit. At 
this point, one can ask two questions: 1) would it be possible 
to build ZAMS models with higher angular momentum content, 
which would remain near the critical limit all along during the 
MS phase? 2) How would such behaviour, if possible, change the 
above results concerning the frequency of stars near the critical 
limit? 

The answer to the first question is yes. It would suffice to dis- 
tribute differently the angular momentum or the angular velocity 
in the ZAMS model, for instance imposing on the ZAMS a non- 
uniform distribution of Q allowing more angular momentum to 
be contained in the inner parts of the star. But this would not 
change the present results concerning the frequencies of stars 
near the critical limit. Indeed, these frequencies are computed 
on the basis of the observed frequencies of stars with given aver- 
aged rotational velocities during the MS phase. Stars with higher 
angular momentum than the most rapid rotators computed here 
will present higher values of the averaged rotational velocities 
during the MS and are thus already implicitly accounted for. 

7.2. Sensitivity on the velocity distribution 

Let us first note that the measure of the rotational velocity 
through the Doppler effect is not free from spurious effects. First 
this technique gives access only to vsini where i is the angle 
between the line of sight and the rotational axis. This quantity 
cannot be known except if some other observations as e.g. a ro- 
tational period obtained from variability induced by the presence 
of dark spot at the surface are available. But even in that case, 
the knowledge of v requires the knowledge of the radius (more 
precisely of the equatorial radius). Moreover if the surface ro- 
tates differentially like the solar surface, the period thus obtained 
might only be partially representative of the surface rotation. 
Second due to the von Zeipel theorem, the fast rotating regions 
of the star, i.e. the zones near the equator, are less luminous than 
the polar regions which have small linear velocities. This effect 
produces an underestimate of the rotational velocities. Townsend 
et al. (2004) obtain the following behaviour between the FWHM 
of the Hel 4471 line for a B2-type star and the equatorial ve- 



locity: the relation is nearly perfectly linear for sin i — 1 and 
v/v C rit < 0.7 (for a standard metallicity 9 Mq star, this cor- 
responds to velocities inferior to 300 km s _1 on the ZAMS). 
In that velocity regime, the FWHM measurements, due to the 
above effects, underestimate the velocities by at most 3-4%. For 
v/Vcat > 0.7, the relation is no longer linear, the FWHM in- 
creases much more slowly with the velocity, it even goes through 
a maximum, decreasing a little for velocities near the critical 
limit. In that case, the FWHM becomes a poor indicator of the 
surface velocity, a star rotating at the critical limit presenting the 
same line broadening as a star rotating at 80% the critical veloc- 
ity. Said in other words, this technique may underestimate the 
real velocity of the star by 20-25% for very fast rotators (typi- 
cally for v/v a i t > 0.75, i.e. for velocities superior to about 340 
km s _1 , for a 9Mq at standard metallicity). To study this effect, 
we have done the following experiment: we have modified the 
velocity distribution given by Huang & Gies (2006a) assuming 
that between 300 and 420 km s _1 all the velocities v should be 
enhanced by a factor equal to 1+(l>-300)/120*0.3 {i.e. a linear 
enhancement amounting at most to a factor 1.3 at 420 km/s). 
Using such a modified velocity distribution, shown by the dotted 
line in Fig. 20, we obtain that the fraction of stars at the critical 
limit for an age equal to 25 Myr at Z = 0.020 would be equal 
to 0.055 instead of 0.049 with the original velocity distribution 
of Huang & Gies 2006a (continuous line in Fig. 20). Thus this 
effect would not significantly change the results. 

In order to obtain a fraction of stars at the critical limit of 
about 10% at Z = 0.020, it is necessary to considerably in- 
crease the fraction of stars born with Q/f2 cr i t superior to 0.6 on 
the ZAMS. Typically passing from the original velocity distribu- 
tion (continuous line in Fig. 20) to a uniform distribution (dashed 
line), one obtains that the fraction of stars with values of Q/Q cr i t 
on the ZAMS superior to 0.6, 0.7, 0.8 and 0.9, passes from 37.5, 
26.5, 15.5 and 6.5% to 40, 30, 20 and 10% respectively. Using 
such a uniform velocity distribution, one obtains that the fraction 
of stars at the critical limit equal 12% at Z = 0.020 and age equal 
to 25 Myr (instead of 5% with the original velocity distribution). 

7.3. Sensitivity on mass loss and magnetic fields 

Mass loss, as we have already discussed above, plays a key role 
in preventing the high mass stars to reach the critical limit. We 
showed that using two different laws for the mass loss rate, very 
different behaviours are obtained for the 60 Mq models at stan- 
dard metallicity. It underlines the importance of this ingredient in 
shaping the evolution of the surface velocity in this mass range. 
For lower initial mass stars, typically for 20 Mq, mass loss be- 
comes less and less important (at least in absence of any surface 
magnetic field). 

We can wonder whether mass loss anisotropics would have 
an impact on the results obtained in this work (which assumed 
spherical winds). Let us recall that when a star is rotating at 
f /fcrit - 0.9, winds of hot stars will become much more in- 
tense at the pole than at the equator. This effect would accelerate 
the evolution towards the critical limit. However, accounting for 
this effect would hardly change the evolution we computed here. 
Indeed, v/v cr [ t becomes superior to about 0.8-0.9 only when the 
surface velocity increases anyway in an already very rapid way 
(see Fig. 15). 

A magnetic field, following the Tayler-Spruit dynamo 
(Spruit 2002), will tend to impose solid body rotation during the 
MS phase (Maeder & Meynet 2005), thus making the approach 
of the critical limit easier (stronger coupling between the core 
and the envelope). We can expect that the inclusion of this effect 
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Fig. 22. The empty symbols correspond to the observed fraction of Be 
stars among B0-B3 type objects in clusters with ages ranging from 10 
to 25 Myr. Data are taken from Wisniewski & Bjorkman 2006 (circles), 
Keller et al. 1999 (triangles) and Maeder et al. 1999 (squares). The filled 
symbols correspond to predictions of the given models for Z = 0.002 
and Z = 0.020 and for ages equal to 10 and 25 Myr respectively (we es- 
timated the fraction of Be in a two magnitude interval below the turn off, 
as did Maeder et al. 1999). The hatched zone labelled by I correspond 
to stars predicted at the critical limit. The hatched zone labelled by II 
shows the fractions of stars with u/u a n > 0.7. These two cases were ob- 
tained assuming a velocity distribution of Huang & Gies (2006a). The 
third hatched region, labelled III corresponds to stars with vlv„- A > 0.7 
and assuming that at low metallicity the initial distribution of rotation 
contains faster rotators (see text). 

would have an important impact on the present results, making 
the fraction of stars at the critical limit more important than in 
the present work. It has however to be noted that at present the 
Tayler-Spruit dynamo is still much debated and some more work 
is needed before assessing whether this mechanism can be active 
in stars or not. 

8. Possible links with the Be stars 

There are good indications that the Be star phenomenon is re- 
lated to the fast rotation of the star. For instance Martayan et al. 
(2006a) showed that the initial velocities of the Be stars are sig- 
nificantly higher than the initial velocities of the normal B stars, 
giving some support to the view that only stars with a sufficiently 
high initial velocity can go through a Be episode. Let us note that 
when material is ejected into a Keplerian disk the star loses an- 
gular momentum, slows down and evolves away from the limit. 
Secular evolution will then bring back the star towards this limit. 
Probably the phases in between the loss of material are much 
longer than the phases during which the star loses material into 
the Keplerian disk. Thus most of the Be stars might be observed 
during these quiet phases when they are evolving back towards 
the limit. It means that the rotation of observed Be stars might be 
somewhat lower that the rotation at which outbursts occur. If the 
link between Be stars and fast rotation is now well accepted, it 
is however not clear what rotation rate is needed for a star to be- 
come a Be star, or why this process is restricted mainly to early 
B type stars. Is rotation the only key parameter governing the 
transformation of a star into a Be star? When does such a trans- 
formation occur? In the present section we would like to provide 
some answers to the above questions. 

Comparisons between number fractions of fast rotating stars 
as predicted by the present models (filled symbols and hatched 



regions) and as observed in clusters (empty symbols, see cap- 
tion) is shown in Fig. 22. We see first that if we consider only 
stars that are strictly at the critical limit (hatched zone I), the 
prediction falls well below the observed points, all the more so 
as there are also rapidly rotating stars that are not Be stars, or Be 
stars that are observed between outbursts. Thus according to the 
present stellar models, the number fraction of stars at the critical 
limit is not sufficient to account for the observed fraction of Be 
stars. 

Let us mention that there are two proposed channels for the 
formation of Be stars: 1) the single star channel, where the crit- 
ical limit is reached because of the natural evolution of a sin- 
gle rotating star; and 2) the binarity channel, with the spin-up 
of the member of a binary system through mass transfer from 
its companion. According to McSwain & Gies (2005a,b), about 
75% of the Be stars detected may have been spun-up by bi- 
nary mass transfer, while most of the remaining Be stars were 
likely rapid rotators at birth. There are indeed some known bi- 
naries among the observed Be stars, but in most cases, no bi- 
narity has been yet clearly assessed. It may however be very 
difficult to detect, because of the long orbital period and the 
very small orbital motions, and in some cases, decades of care- 
ful measurements have been necessary before detecting a com- 
panion (see Gies et al. 1998 for (f> Per, Maintz et al. 2005 for 59 
Cyg, Rivinius et al. 2004 for FY CMa). According to Negueruela 
(2007), the binarity channel may be significant, but it is unlikely 
that most Be stars may form through this channel. From a theo- 
retical point of view, Pols et al. ( 1 99 1 ) deduce that probably 40% 
(at most 60%) of the Be stars may form through mass transfer 
spin-up. Consequently, the binary channel alone could account 
for a (Be/(B+Be)) ratio of the order of 8-10% if we take the 
(Be/(B+Be)) value found by Martayan et al. (2006b). In that last 
case, our prediction remains thus too low even when discounting 
the binarity channel contribution. 

In the frame of the single star scenario, if we assume that the 
Be phenomenon may occur already at v/u cr i t > 0.7 (see hatched 
zone II), then the situation improves somewhat, however the the- 
oretical results do not show any increase of the fraction of Be 
stars with decreasing metallicity, as some observations indicate 
(Maeder et al. 1999; Wisniewski & Bjorkman 2006). It might 
be that part of this trend is due to an observational bias. For in- 
stance McSwain & Gies (2005b) found that Be stars were gener- 
ally found among the brighter B-stars in their survey of southern 
sky Galactic clusters. A similar trend was found by Keller et al. 
(2001) in h andx Per. Since metal poor clusters are more distant, 
observations will pick up mainly the brighter end of the B stars 
sample where the fraction of Be stars is higher. However, it is 
not clear that such a bias will have a great impact when clusters 
of similar ages are compared: provided the observed samples are 
complete and the count is performed in the same absolute magni- 
tude interval, then the observed number ratios at different metal- 
licities can be compared. To reproduce the observed trend with 
the metallicity obtained by Maeder et al. (1999) and more re- 
cently by Wisniewski & Bjorkman (2006) with the present mod- 
els, we have to assume that the fraction of fast rotators at birth is 
higher in metal poor regions than in metal rich one. As a numer- 
ical example, supposing that the initial distribution of rotation 
at Z = 0.002 is given by a uniform distribution (dashed line in 
Fig. 20), we then obtain the hatched zone III, which agrees rea- 
sonably well with the observations. From the discussion above, 
we deduce, that to reproduce the observational number fraction 
of Be stars, first we have to relax the hypothesis that Be stars 
are strictly at the critical limit: the Be phenomenon should occur 
already for v/v CT i t ^ 0.7; second, the distribution of initial val- 
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ues of Q/Q cr it on the ZAMS at low metallicity should contain a 
higher fraction of fast rotating stars. 

Let us discuss the two above conclusions. The first conclu- 
sion would tell us that, although fast rotation plays an important 
role in transforming the star into a Be star, it is not sufficient 
alone. Above some value of u/v CI it some other processes, like 
for instance pulsations, must help in launching material into a 
disk. 

Now what about the second conclusion: a higher fraction of 
fast rotators at low Z? Keller (2004) measured the rotational ve- 
locity of 100 MS early B-type stars in clusters of the LMC and 
compared the results with observations of early B-type stars in 
clusters of the solar neighbourhood. He obtains that the LMC 
stars are faster rotators than the galactic ones: the mean values 
of usin; is 116 km s _1 for the galactic stars and 146 km s _1 
for the LMC stars. Part of the difference can be due to the dif- 
ference in mass loss between galactic and LMC stars, but part 
might also be due to differences in the initial distribution of the 
velocities. Penny et al. (2004), on the other hand, find no dif- 
ference between the velocities of O-type stars in the Magellanic 
Clouds and in the Galaxy. Thus the picture is at present not clear 
and probably new surveys, as e.g. the VLT-Flames survey will 
provide very soon new insights on that question. At the moment 
we still consider the question as opened. Present results would 
however clearly favour a higher fraction of stars born with high 
values of £2/Q C rit at low metallicities. 

A point which the models have to explain is when the Be 
phenomenon does appear during the MS phase. Unfortunately 
no clear picture emerges from the observations: some authors 
find the highest fraction of Be stars among the spectral types B0- 
B2, appearing most often at ages between 13-25 Myr (Fabregat 
& Torrejon 2000). Recently Fremat et al. (2006) who determined 
the fundamental parameters of the Be stars which would be ob- 
servable by COROT found that the main part of their stellar 
sample falls in the second half of the MS lifetime. This would 
favour a view according to which the Be star phenomenon ap- 
pears preferentially at the end of the MS phase. Other authors 
(as e.g. Zorec & Briot 1997) conclude that the Be phenomenon 
is not associated with a particular evolutionary phase and can 
occur during the whole MS phase, at least for some range of 
initial masses (see Zorec et al. 2005). With the adopted distribu- 
tion of initial velocities, present models predict that most near 
critical velocity stars will be found at evolved stages during the 
MS phase and would thus provide an explanation for the obser- 
vation of Fabregat & Torrejon (2000). From a physical point of 
view, we think that the Be star phenomenon may be limited to a 
given spectral type range for the following reasons: for the ear- 
lier spectral types than the favoured one, mass loss prevents the 
star from approaching the critical limit and for the later spec- 
tral types, the meridional currents are not efficient enough for 
accelerating the envelope. Models would predict a shift towards 
earlier spectral type when the metallicity decreases due to the 
decrease of mass loss with the metallicity (extending the range 
towards earlier spectral type) and the decrease of the meridional 
velocity in denser layers (increasing the minimum initial mass 
for meridional currents being efficient enough to accelerate the 
surface). Of course, beyond these physical reasons are the effects 
of a possible variation of the distribution of initial rotation as a 
function of metallicity, which would somewhat blur the above 
picture. And now if we relax the condition for becoming a Be 
star to v/ u cr [ t ^ 0.7, we allow this phenomenon to appear earlier 
in the life of a star. 

The present models predict the approach of high values of 
v/Vcm (0.8-0.99) only when the star is already evolved into the 



MS phase. As mentioned above, the use of a different internal 
profile of Q in our ZAMS models would allow us to explore still 
more extreme cases of rotation, where the star may reach the crit- 
ical limit, or at least remain near to it, from the beginning of the 
MS phase. However these stars should be very rare. The neces- 
sity for the stars to be somewhat evolved in order to reach high 
surface velocities implies that these stars should present at their 
surface N-enrichments. These surface enrichments might how- 
ever be very small for stars in the lower part of the mass range 
possibly evolving into a Be star phase (see the cases of our 3 Mq 
stellar model at Z = 0.020 in Fig. 17). Observations of N surface 
enrichments for Be stars are still scarce, probably because of the 
difficulty of determining abundances from greatly widened ab- 
sorption lines. Fremat et al. (2004) obtain N/C ratios for four Be 
stars which are clearly higher than what is observed at the sur- 
face of slowly rotating B-type stars. They obtain (N/C)/(N/C)i„i 
of the order of a factor 2, which would well correspond to what 
is expected for near critically rotating 9 Mq stars at at Z = 0.020 
(see middle panel of Fig. 17). 

9. Conclusions 

From the above computations we can deduce that 

1 . As long as the mass loss rates are not too important, merid- 
ional currents, by advecting angular momentum from the in- 
ner regions to the outer layers, bring the v/v CI i t ratio close to 
1 during the MS phase. 

2. Since the velocity of the meridional currents in the outer lay- 
ers scales with the inverse of the density, the process be- 
comes more efficient for stars of higher initial mass and/or 
higher initial metallicity. 

3. At high metallicity however, mass loss becomes more and 
more important and can prevent the stars from reaching the 
break-up limit. 

4. In stellar clusters, one expects that the fraction of stars hav- 
ing v/v CI it ^ 0.8 becomes maximum for ages between 20-32 
Myr at Z = 0.020. This range of ages shifts to older ages at 
higher metallicities and to younger ages at lower metallici- 
ties. Such a shift results from the dependence on the metal- 
licity and the mass of the intensity of the mass loss rates and 
of the velocities of the meridional currents. If we relax the 
condition and look at stars with i//f C rit ~ 0.7, the maximum 
number will be found in clusters of ages around 10 Myr and 
below, and this at all non-zero metallicities. 

5. Be stars might be the natural outcome of stars with an ini- 
tial rotational velocity in the upper tail of the initial velocity 
distribution. To explain their frequency from the present ro- 
tating stellar models, minimum values of v/v cr [ t <: 0.7 have 
to be reached during the MS for the Be star phenomenon to 
occur. 

6. Depending on when this limit is reached, one expects more 
or less high surface enrichments. If the limit is reached very 
early during the MS phase, no enrichment is expected, while 
if this limit is reached at the end of the MS phase, high N/C 
and N/O ratios are expected. For Be stars originating from 
small initial mass stars (typically from 3 Mq at standard 
metallicity) one expects no or small surface enrichments. 

7. To reproduce the higher fractions of Be stars at low metal- 
licity, we must assume that in metal poor regions, a higher 
number of stars are born with high values of Q/£2 cr i t . 

In a future work, it will be interesting to explore the effects of a 
magnetic field. As already mentioned above, this will probably 
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change significantly the results presented here, allowing the stars 
to approach more easily the critical limit. It will be interesting to 
see if, in that case, it would still be necessary to invoke a higher 
proportion of fast rotators at low Z. 
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